Using H(z) data as a probe of the concordance model 
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Direct observations of the Hubble rate, from cosmic chronometers and the radial baryon acoustic 
oscillation scale, can out-perform supernovae observations in understanding the expansion history, 
because supernovae observations need to be differentiated to extract H(z). We use existing H(z) 
data and smooth the data using a new Gaussian Processes package, GaPP, from which we can 
also estimate derivatives. The obtained Hubble rate and its derivatives are used to reconstruct the 
equation of state of dark energy and to perform consistency tests of the ACDM model, some of 
which are newly devised here. Current data are consistent with the concordance model, but are 
rather sparse. Future observations will provide a dramatic improvement in our ability to constrain 
or refute the concordance model of cosmology. We produce simulated data to illustrate how effective 
H(z) data will be in combination with Gaussian Processes. 



I. INTRODUCTION 

Different methods and data sets are being used to re- 
construct the dark energy (DE) equation of state w = 
Pde/ P&e and thereby also to test the concordance model 
(which has w = —1). The results vary significantly ac- 
cording to the methods and data sets used, and the error 
bars and uncertainties are large. It is clear that higher- 
precision data are needed for an effective reconstruction 
and for robust testing of models. But just as important, 
more effort is needed to improve the statistical meth- 
ods and the design of observational tests. In particular, 
there is a need for effective model-independent statistical 
methods and for tests that target the concordance model. 

One of the most direct ways to reconstruct w is via 
supernovae (SNIa) observations that give the luminosity 
distance c?l- Model- independent approaches to recon- 
structing w have been developed [1H19|. SNIa observa- 
tions lead indirectly to H(z) via the derivative d' L (z). 
Then we need the second derivative of dz,(z) to recon- 
struct w. This is very challenging for any reconstruction 
technique since any noise on the measured dz,(z) will be 
magnified in the derivatives. The problem can be less- 
ened if direct H{z) data are used because only the first 
derivative needs to be calculated to determine w{z). 

In this paper we focus on observations that directly 
give H(z). Presently, this may be derived from differ- 
ential ages of galaxies ('cosmic chronometers') and from 
the radial baryon acoustic oscillation (BAO) scale in the 
galaxy distribution. Compared to SNIa observations, less 
H(z) observational data are needed to reconstruct w with 
the same accuracy. For the cosmic chronometer data, it 
has been estimated [2(| that 64 data points with the ac- 
curacy of the measurements in [2l| are needed to achieve 
the same reconstruction accuracy as from the Constitu- 
tion SNIa data (22| . 

We use a model-independent method for smoothing 
H(z) data to also perform consistency tests of the con- 
cordance model (flat ACDM) and of curved ACDM mod- 



els. These consistency tests are formulated as functions 
of H(z) and its derivatives which are constant or zero 
in ACDM, independently of the parameters of the model 
(see [23[ for a review). Deviations from a constant func- 
tion indicate problems with our assumptions about dark 
energy, theory of gravity, or perhaps something else, but 
without the usual problems of postulating an alternative 
to ACDM. Some of the tests we use here are given for 
the first time. 

Gaussian processes (CP) provide a model-independent 
smoothing technique that can meet the challenges of re- 
constructing derivatives from data [24|,[25|]. We follow the 
same GP approach that has been applied to supernova 
data in a previous work [n| by some of the authors of this 
paper. We use IGaPPl (Gaussian Processes in Python), 
their publicly available code 1 . (See [3, EH for different 
uses of GP in this context.) A brief description of the 
GP algorithm is given in Appendix [A"l 



II. TESTING ACDM 



The Friedmann equation, 



H 2 (z) 



Hi 



(1 — Q m — VLk) exp 



n m (i + zf + n K {i + z) 2 

1 + 



l + z' 



-dz' 



can be rearranged to give 



w{z) 



Pde _ 2(1 + z)hti - 3h 2 + Cl K (l + zf 

pde ~ 3[h 2 -n m (i + z ) 3 -n K (i + zy 



,(i) 



(2) 



In principle, given h(z) data we can smooth it, attempt to 
estimate its derivative, and reconstruct w(z). However, 



1 http: //www. aege .uct . ac . za/~ seikel/GAPP/index.html 
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FIG. 1: h(z) — H(z)/Ho (top) and h'(z) (bottom) reconstructed from cosmic chronometer data (left), BAO data (middle) and 
CC+BAO data (right), using Gaussian processes. Shaded areas represent 68% and 95% confidence levels. The dashed (red) 
curve is flat ACDM with Q m = 0.27; the solid (blue) curve is the GP mean. Note that while the BAO data appear to give an 
inconsistent h'(z), this is driven by the two highest redshift points both of which happen to lie below the flat ACDM curve. 



reconstruction of w(z) is compromised by various difficul- 
ties. It depends on the values of fi m and Qk , so we need 
independent information about these parameters when 
we reconstruct w{z) from H(z) data. These are difficult 
to estimate without assuming a form for w(z) [26l - [28j . 

These difficulties reflect the fact that we cannot use 
data to construct physical models - rather, we need to use 
data to test physical models. The ACDM model could 
be tested by looking for deviations from w = — 1. How- 
ever, there is a more focused approach: to develop null 
hypotheses for ACDM, independently of the parameters 
Cl m and VLk HH. 

To test the concordance model - i.e. flat ACDM - 
we can use jf]) to define a diagnostic function of redshift 

MM- 



h 2 -i 



z(3 + 3z + z 2 )' 



(3) 



Then 



0^' (z) = Ct m implies the concordance model. 

If Om\z) is not a constant, this is a signal of an al- 
ternative dark energy or modified gravity model. Given 
observed h(z) data, we can estimate confidence limits for 
Om • If these are not consistent with a constant value, 
we can rule out the concordance model. 

It is more effective to measure deviations from zero 
than from a constant. The more effective diagnostic is 
thus the vanishing of the derivative Om (z). This is 
equivalent to = 0, where [29| 

C {1) = 3(1 + zf{l- h 2 ) + 2z(3 + 3z + z 2 )hti. (4) 



The null test is therefore 

O- 1 ' ^ falsifies the concordance model. 

To apply this test, we need to reconstruct h'(z) from the 
data. 

If the concordance model is ruled out, it is still possi- 
ble that a curved ACDM model describes the Universe. 
Equations ([1) and (with w = —I) form a linear sys- 
tem for Q rn and VLk- Solving for these parameters we 
can define 

_ a (l + «)(l-tf) + f + (5) 
z z (l + z){3 + z) 

3(1 + z) 2 (h 2 — I) — 2z(3 + 3z + z 2 )hh' 

°K{Z) = T — i r— r ,(6) 



and we have 



0${z) 
O k (z) 



z 2 (l + z)(3 + z) 

fi m implies ACDM, 
Q K implies ACDM. 



These quantities are equivalent to those derived in [32| 
in terms of D(z), the dimcnsionlcss comoving luminosity 
distance. The D(z) forms contain second derivatives D" 
whereas the h{z) forms above contain only first deriva- 
tives h! . Given observed Hubble rate data from which we 
can estimate the derivative h'{z), we can then estimate 
confidence limits for O m (z) and O k (z). If these are not 
consistent with a constant value, we can rule out ACDM 
in general, and conclude that dark energy has w 1 (or 
there is modified gravity). 

The more effective diagnostic of these consistency tests 
is the vanishing of the derivatives of (O and ((5]). The 
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vanishing of Om is equivalent to C}- 2 ^ = 0, where 

£ (2) {z) = 3(l + z) 2 (h 2 -l)-2z(3 + 6z + 2z 2 )hti 

+ z 2 (3 + z){l + z)(h ,2 + hh"). (7) 

Then 

£ (2) (z) ^ falsifies ACDM. 

The vanishing of O K does not give any independent 
information - it is also equivalent to £' 2 ' = 0. 

Given observations of h(z), we can construct this func- 
tion independently of the parameters of the model and 
test ACDM by measuring consistency with zero. This 
has the advantage that it is easier to detect deviations 
from zero rather than a constant, but at the expense 
of requiring an extra derivative in the observable. This 
is akin to detecting deviations from constant in w, but 
without reliance on the parameters of the model. 

For the application of these consistency tests, it is cru- 
cial to use a model-independent method to reconstruct 
o£\ {2 \ O k , and C^l Model-dependent ap- 

proaches have the problem that they affect or even deter- 
mine the outcome of the consistency test: While fitting 
a ACDM model to the data would always lead to a re- 
sult that is consistent with ACDM, fitting a model that 
does not include ACDM as a special case would result in 
inconsistencies with ACDM. The only mo del- dependent 
approches that do not entirely determine the outcome 
of the test are those assuming a model which includes 
ACDM as a special case. Nevertheless, they affect the 
result by forcing the data into a specific parametrisation, 
which might not reflect the true model. The only way to 
avoid this problem is to use a non-parametric approach. 
Here, we use Gaussian processes, which are described in 
Appendix 

III. RECONSTRUCTION AND CONSISTENCY 
TESTS FROM H(z) DATA 



The second way to measure H(z) is the observed line- 
of-sight redshift separation Az of the baryonic acoustic 
oscillation (BAO) feature in the galaxy 2-point correla- 
tion function [37HI1], 



where r s (zd) is the sound horizon at the baryon drag 
epoch. 



Results: real data 

We use the following H(z) data sets: 

CC: 18 cosmic chronometer data points [3a |. 

BAO: 6 radial BAO data points 

CC+BAO: Combination of CC and BAO sets. 

We normalize H(z) using Ho = 70.4 ± 
2.5kms~ 1 Mpc~ 1 . The uncertainty in Ho is trans- 
ferred to h(z) as a 2 = (a 2 H /H 2 ) + (H 2 /H^)a 2 Ho . The 
reconstructed functions h(z) and h'(z) are shown in 
Fig. [TJ The shaded regions correspond to the 68% and 
95% confidence levels (CL). The true model is expected 
to lie 68% of the plotted redshift range within the 68% 
CL. Note that this is only an expectation value. The 
actual value for a specific function may deviate from the 
expectation. The dependence of the actual percentage 
on the smoothness of the function has been analysed in 

Figure [2] shows the reconstruction of Om • The re- 
construction of O m and Ok is shown in Fig. [21 and 
Fig. [3] gives C^ 1 ' and O- 2 '. We actually plot a mod- 
ified C m = C/(l + z) 6 which stabilises the errors at 
high redshift without affecting the consistency condi- 
tion. The reconstructed w(z), also requiring h' , is shown 
in Fig. 21 where we assume the concordance values 
O m = 0.275 ± 0.016 and J] K = O0. 



Cosmic chronometers are based on observations of the 
differential ages of galaxies [2l|, [33T - I351 ] . The Hubble rate 
at an emitter with redshift z is 



H(z) = 



1 dz 

1 + z dtp 



(8) 



where t e is the proper time of emission. The differential 
method uses passively evolving galaxies formed at the 
same time to determine the age difference At e in a small 
redshift bin Az, assuming a Fricdmann background. To 
find old galaxies sharing the same formation time, we 
have to look for the oldest stars in both galaxies and show 
that they have the same age. This method is effective; 
but while the differential approach significantly reduces 
the systematics that would be present when determining 
the absolute ages of galaxies, it still faces uncertainties 
due to the assumptions that are made to estimate the 
age. 



Results: mock data 

To demonstrate how a larger number of data will affect 
our results when reconstructing w and testing ACDM, 
we simulated a data set of 64 points for H(z), drawing 
the error from a Gaussian distribution Af(a, e) with a = 
10.64z + 8.86 and e = 0.125(12.46z + 3.23), adapting the 
methodology of [2(| ■ 

We simulated data points for two different models: 
Concordance model, £Ik = 0, Q m = 0.27. 
A model with slowly evolving equation of state: 



11 / 1 

w{z) = -- + -tanh3^z - - 



and the same concordance density parameters. 
The GP reconstructions arc shown in Figs. [SH5J 



(10) 
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FIG. 2: Om (z) (top), (z) (middle) and Ok(z) (bottom) reconstructed from cosmic chronometers (left), BAO (middle) 
and CC+BAO (right). For o£\z), the dashed (red) curve is flat ACDM. For 0$ (z) and Ok(z) it is a curved ACDM model. 
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FIG. 3: 
(middle) and 



= £ (1) /(1 + zf (top) and d$ = £ {2) /(l + zf (bottom) reconstructed from cosmic chronometers (left), BAO 
CC+BAO (right). The dashed (red) curve is a ACDM model. 
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FIG. 5: h(z) (top), h (z) (middle) and h"(z) (bottom) reconstructed from simulated data, assuming a concordance model (left) 
and model (|10[) with slowly evolving w(z) (right). 



Discussion 

Figure^ shows that for the CC and CC+BAO data (18 
and 24 points), we get good reconstructions when there 
is no differentiation of h(z) involved. The BAO data set 
only contains 6 data points up to redshift 0.73. Beyond 
that redshift, the reconstruction differs significantly from 
ACDM. The results from the CC and CC+BAO sets are 
however in very good agreement with ACDM. 

The BAO data appear to be inconsistent with the con- 
cordance model. However, 6 data points are not suffi- 



cient for a reliable reconstruction. The two data points 
with highest redshift happen to be below the concordance 
curve, which pulls the reconstructed curve down. This 
is probably just a coincidence, but it illustrates the im- 
portance of having the derivative of the data consistent 
with the model, as well as the data itself. Current and 
upcoming large- volume surveys, such as BOSS |4l|, EU- 
CLID g2( and SKA H|, will provide radial BAO mea- 
surements of increasing number and precision. 

The reconstruction of Om and Ok shown in Fig. [5] 
is more challenging for the available data set, since we 
need the first derivative of h. With present data sets, the 
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FIG. 6: Om\z) (top), Om\z) (middle) and Ok{z) (bottom) reconstructed from simulated data, assuming a concordance model 
(left) and model (right). 
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uncertainties in the reconstruction are quite large. Using 
CC and CC+BAO, these results as well as the results for 
£W and shown in Fig. 02 are consistent with ACDM. 

For the mock data sets, Figs. [5] and [6] show that the GP 
reconstructions recovers the assumed models very effec- 
tively. We can clearly distinguish the model with slowly 
evolving w(z) from ACDM in 0$ . For 0$ and K , 
the reconstruction errors are too large to see this differ- 
ence. The same is true for consistency tests and 
shown in Fig. [7J 

The reconstruction of the equation of state w(z) also 
shows a clear difference of the two models, assuming we 
can accurately determine Hq, fl m and fix separately 
from w(z): see Fig. [HI GP works very well to recover the 
assumed if. With less than 100 data points, we can re- 
construct a dynamical dark energy model far better than 
is achievable using thousands of SNIa data - compare to 
analogous reconstructions in 



IV. CONCLUSIONS 

We have considered the information that current and 
future H(z) data can give us. Currently such data come 
from cosmic chronometers and BAO data, and is plainly 
consistent with the concordance model. Future data, 
however, will provide a powerful discriminator between 
models. It is remarkable how few data points are required 
compared to supcrnovac: to reconstruct w(z) accurately 
in our non-parametric way requires many thousands of 
SNIa, compared to less than 100 H(z) data points. 

We have derived and analysed new consistency tests for 
the ACDM model, which we have formulated in terms of 
H(z) directly, rather than using the more familiar dis- 
tance function [23|, [HJ . By smoothing the data points 
using Gaussian process, we have shown that these can be 
very effective in determining that ACDM is the incorrect 
model, but without having to assume the key parameters 
VL m and which currently only have constraints de- 
rived by assuming ACDM or a similar alternative. These 
tests not only require that the data points themselves 



are consistent with the model, but that their derivative 
is also. 

Future data which directly measures the expansion his- 
tory will therefore play an important role in future dark 
energy studies. 
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Appendix A: Gaussian Processes 

For a data set {{zi,yi)\i = l,...,n}, where Z rep- 
resents the training points z,-, i.e. the locations of the 
observations, we want to reconstruct the function that 
describes the data at the test input points Z* . 

A Gaussian Process is a distribution over functions and 
is thus a generalization of a Gaussian distribution. It is 
defined by the mean n(z) and covariance k(z,z): 



f(z)~gV(fi(z),k(z,z)) 



(Al) 



At each Zj, the value /(zj) is drawn from a Gaus- 
sian distribution with mean /i(z,) and variance k(zi, Z,-). 
/(zj) and f(zj) are correlated by the covariance function 

fe(Z}, Zj). 

Choosing the covariance function is one of the main 
points for achieving satisfactory results. The squared ex- 
ponential is a general purpose covariance function, which 
we use throughout this paper: 



k(zi, Zj) = aj exp 



2£ 2 



(A2) 



The 'hyperparameters' are <7/ (signal variance) and I 
(characteristic length scale), t can be thought of as the 
distance moved in input space before the function value 



8 



changes significantly. Of describes the typical change in 
y-direction. In contrast to actual parameters, they do 
not specify the exact form of a function, but describe 
typical changes in the function value. 

For Z* i the covariance matrix is given by 
[K(Z*,Z*)]ij = k(z*,z*). Then the vector /* with 
entries f(z*) is drawn from a Gaussian distribution: 



/* ~JV(n{Z*),K{Z*,Z*)). 



(A3) 



This can be considered as a prior for the distribution 
of /*. One needs to add observational information to 
obtain the posterior distribution. 

The observational data have a covariance matrix C. 
For uncorrelated data, C is a diagonal matrix with en- 
tries Ui. The combined distribution for /* and the ob- 
servations y is given by: 



V 

r 



K(Z,Z) + C K(Z,Z* 
K(Z*,Z) K(Z*,Z* 



(A4) 



While the values of y are already known, we want to 
reconstruct /* . Thus, we are interested in the conditional 
distribution 



(A5) 



where 



/* = n* + K(Z*, Z) [K(Z, Z) + C\-\y - M ) (A6) 
cov(/*) = K{Z*,Z*) 

- K(Z*, Z) [K(Z, Z) + C}- 1 K(Z, Z*), (A7) 

are the mean and covariance of /*, respectively. The 
variance of /* is simply the diagonal of cov(/*). Equa- 
tion (|A5[) is the posterior distribution of the function 
given the data and the prior (|A3|) . 

In order to use this equation, we need to know the 
values of the hyperparameters <jf and I. They can be 
trained by maximizing the log marginal likelihood: 



ln£ = \np(y\Z 1 af,t) 

= - l -{y-v) T [K{Z,Z) 



CY l {y-n) 



I 77 

- In \K(Z, Z) + C\ - - In 2tt . (A8) 



Note that this likelihood only depends on the observa- 
tional data, but is independent of the locations Z* where 
the function is to be reconstructed. 

Derivatives of the function can be reconstructed in a 
similar way. For the first derivative, the conditional dis- 
tribution is given by [19 ]: 



f*'\Z*,Z,y~M(t', C ov(f*')) , 



(A9) 



where 



/*' = M *' + K'(Z*,Z) [K(Z, Z) + C]-\y - M ) (A10) 
cov(/*') = K"(Z*,Z*) 

- K'(Z*, Z) [K{Z, Z) + C]" 1 K'(Z, Z*). (All) 
For the covariance matrices, we use the notation: 



[K'(Z,Z*)h = 
\K"{Z\Z*)\ l3 = 



dk{z l ,z* j 



d 2 Kzt,z*) 
dz* dz* 



(A12) 
(A13) 



K'(Z*,Z) is the transpose of K'(Z, Z*). 

To calculate a function g(f, f) which depends on / 
and /', we also need to know the covariances between 
/* = f{z*) and /*' = f'(z*) at each point z* where g is 
to be reconstructed. This covariance is given by: 



ewer,/*') 



dk(z*,S) 



dz 



(A14) 



K'(z*,Z)[K(Z,Z) + C] 1 K{Z,z*). 



g* = g{z*) is then determined by Monte Carlo sampling, 
where in each step /* and /*' are drawn from a multi- 
variate normal distribution: 



/* 
/*' 



' Af 



f* 
/*' 



var(/*) cov(/*,/*') 
cov(/*,/*') var(/*') 



(A15) 
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